Single light scattering: Bubbles versus droplets 
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The intensity and polarization characteristics of scattered light beams are studied for bubbles in 
water and compared to the case of water droplets in air. The size distributions of the scatters in both 
systems are assumed to be the same. It is shown that spherical polydispersions of water droplets give 
generally larger depolarization effects than bubbles, leading to larger entropy production by droplets 


than bubbles. 
[DOI: 10.1119/1.1621030] 


I. INTRODUCTION 


The interaction of light with macroscopic liquid or solid 
particles can be studied in the framework of Mie theory.! 
Usually one is interested in the scattered light intensity and 
polarization characteristics as functions of the scattering 
angle, particle sizes or shapes, and their refractive indices.”~ 
The goal of this paper is to identify differences in light scat- 
tering that are due to different types of mixtures of the same 
substances. 

In particular, we consider the case of air bubbles in water 
and water droplets in air. Our results are of interest for stud- 
ies of light propagation through a cloudy atmosphere and the 
ocean. In particular, bubble clouds play an essential role in 
radiative transport close to the ocean surface, especially if a 
strong wind is present. Above the ocean surface there is a 
cloud of droplets and below it bubble clouds develop. 

Droplets in clouds and bubbles in water are usually much 
larger than the wavelength of visible light. They also are 
polydisperse. Our calculations are for spherical polydisper- 
sions with the following particle size distribution, which 
is frequently used for studies of light propagation in 
clouds:* 


f(a) =Aa? exp(—0.6a), (1) 


where A =35/(24X 58), which insures that 


| faraa=1. (2) 


Asymmetric distributions such as those given by Eq. (1) are 
typical for water droplets in terrestrial clouds.** Note that the 
mode radius ag in Eq. (1) is equal to 10 um. The mode 
radius is defined in such a way that the function in Eq. (1) 
has a maximum at this radius. Note that ag is much larger 
than the wavelength of the incident light à used in the cal- 
culations (A =0.55 um). 


Il. THEORY 


The problem of electromagnetic wave interaction with a 
single macroscopic spherical particle (for example, a cloud 
droplet) was solved roughly a century ago (see, for example, 
Ref. 1). The solution in the far field can be written in terms 
of the Stokes vector parameter of the scattered radiation 5 
(S1:S2:8384), 
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S= ar So, (3) 


where 0 is the scattering angle, k=27/N, N is the wave- 


length, r is the distance to the observation point, and S o is the 
Stokes vector parameter of the incident light. Note that the 
components of the Stokes vector parameter are defined as 
bilinear combinations of the correspondent electric field vec- 
tor E=E,i +E,j.2-4 The vector k=j Xi gives the direction 
of the electromagnetic wave propagation. Namely, we 
have*"* §,;=E,E*+E,Ex, S,=E,E*—E,E% , S,=E,Ex 
+EŤE,, Sg=i(E|E}— EE), where we neglected a com- 
mon multiplier. In electromagnetic scattering problems, it is 
usually assumed that the vector i is in the plane of scattering, 
and the vector j is in the plane perpendicular to the plane of 
scattering. 

The nonzero elements of the 4X4 scattering matrix Ê for 
spheres can be calculated using the following equations:”~ 


P= Py = PIA nlt], 


4 
D1.=P21= A ul? |A 2017], (4) 

D33= O44=Re[A AJ], 

D34= — By3=In(A | A355], ©) 
where A,; and A>, are elements of the amplitude scattering 
matrix 

= 6 
| 0 a k 
and?-* 
” 2j+1 
ay > jor 
“0 (7) 
2j+ 
= + 

An= >, jG+)) la, jm; +b;7;]. 

Here m; is the derivative of the Legendre polynomial 


P (cos 0) and 7; is the derivative of the associated Legendre 
polynomial Pj (cos 0) with respect to cos 6. The amplitude 
coefficients a; and b; are given by complex combinations of 
the spherical Bessel functions, which depend on the complex 
refractive index of particles and the ratio a/ .? Both Leg- 
endre polynomials and spherical Bessel functions can be 
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Table I. The physical meaning of the normalized scattering matrix elements [see Eq. (8)]. 


Stokes vector of the 
incident light 


Stokes vector of the 


scattered light for spherical 


polydispersion 


Degree of polarization 


Unpolarized light 


1 
0 
0 
0 
Right circularly 
polarized light 
1 
0 
0 
1 
Linearly polarized light 
with the azimuth —45° 


Partially linearly 
polarized light 
Py 
Po 
0 
0 


Partially circularly 
polarized light 
Pı; 

0 
0 
Px 


Partially elliptically 
polarized light 


Linear polarization 


Circular polarization 


Dy, 
jae 
Pe Diy 


Circular polarization 


1 Pi » Pa 
0 bı í EO 
=1 —®,, 
0 P34 
found using standard recurrence formulas. The correspond- Ia 
ing equations can be found in Refs. 2—4. Ca a (27+1)(\a;|?+|b,|?). (11) 
j=l 


It should be stressed that the elements of the matrix Ê are 
easily measured using different illumination conditions as 
summarized in Table I. It follows from Table I that the ele- 


ments of the normalized scattering matrix p=®/®,, have a 
clear physical meaning. Namely, it follows that 


1 — Pi 0 (0) 


Pa o e opel 7 
0 0° =p% Pe 


where p; is the degree of linear polarization of scattered light 
assuming the illumination of a particle by unpolarized light 
and p, and p* are the degrees of circular polarization of 
scattered light assuming circularly or linearly (with the azi- 
muth equal to —45°) polarized incident light, respectively.* 
The matrix (8) plays an important role in radiative transfer 
problems.** We will consider the elements of this matrix in 
the following. 

Another important characteristic is the total scattered in- 
tensity for the illumination of a particle by unpolarized light. 
Its angular distribution is given by ®,, [see Eq. (3) and Table 
I]. We will consider, however, not ®,,, but the phase func- 
tion 


47®,,(0) 
P( = ea (9) 
where 
2m (7 
Cu | ®,,(0)sin 0d0 (10) 
0 


is the scattering cross section.’~* The integral (10) can be 
calculated analytically. The result is? 
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Note that it follows by definition: $f7p(0)sin 0d0=1. 
Therefore, the phase function gives the conditional probabil- 
ity of light scattering in a fixed direction.* 

Equations (4)—(7), (9)—(11) are valid for monodisperse 
scatterers only. The conversion to the polydispersed case is 
straightforward.’ We need to substitute ®;; and Csca by the 
average values (®;;) and (C,.,). The angular brackets have 
the following meaning: 


(v)= | Y@ntayaa. (12) 


Ill. RESULTS OF NUMERICAL CALCULATIONS 


The results of calculations of the phase function p(@), 
which gives the conditional probability of light scattering in 
the direction specified by the scattering angle 0, are shown in 
Fig. 1(a). This result was found using Eqs. (9), (4), (11), and 
(7). The function ®,,(@) and the value of Csa in Eq. (9) 
were averaged according to Eq. (12) with the particle size 
distribution (1). 

We see that p( 0) is quantitatively and qualitatively differ- 
ent for bubbles and droplets. In particular there is no rainbow 
scattering for bubbles as there is with droplets (see the angle 
range 120°—160°). There is an increase in the phase function 
as 0— m for droplets, but this increase is not as pronounced 
for bubbles. The probability of photon scattering by bubbles 
in the region from 120° to 180° is considerably smaller for 
bubbles than droplets. Thus, the backscattered signal from 
bubbles is comparatively weak. On the other hand, bubbles 
are stronger scatterers at the lateral angles (from 40° to 
110°). Due to the same size of the droplets and bubbles in 
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Fig. 1. (a) The phase function of water droplets in air (solid line) and air 
bubbles in water (dashed line). The particle size distribution is given by Eq. 
(1). The wavelength is equal to 0.55 um. (b) The phase function of air 
bubbles in water, calculated using the geometric optics approach for a single 
bubble (Ref. 2), with a radius 1 mm (solid line). The results given by the 
dashed lines in (a) and (b) coincide. 


our model, their phase functions coincide when 0< 04, 
where 0,=5° in the case studied. The range [0,04] specifies 
the diffraction zone. 

There is an interesting feature in the behavior of the phase 
function of bubbles when the scattering angle is close to 70°. 
This feature is related to the fact that rays incident on spheri- 
cal particles (with refractive indices n<1) at angles larger 
than the critical angle g=arcsin(n) do not penetrate inside 
the particles.° The critical scattering angle 6, is related to the 
critical angle œ by? 0.= 7-2. This gives for the critical 
scattering angle 0,=82.82° for n=0.75. Light scattering at 
01< 0< 0, is mostly due to the total internal light reflection 
from the surfaces of the bubbles. Detailed studies of wave 
effects at angles close to 0, [see Fig. 1(a)] were performed 
by Fiedler-Ferrari et al.’ 

We compare our calculated results using geometric optics” 
(valid for a/N— œ) for the phase function of bubbles to 
those using wave optics in Fig. l(b). This comparison sug- 
gests that the case studied (ay=10 wm) is very close to the 
asymptotic geometric optics scattering regime in the back- 
ward hemisphere. However, there are differences in the re- 
gion near 0., where total internal reflection takes place.” The 


260 Am. J. Phys., Vol. 72, No. 2, February 2004 


droplets 


bubbles a 


y 


degree of linear polarization 


O 20 40 60 80 100 120 140 160 180 
scattering angle, degrees 


Fig. 2. The angular dependence of the degree of linear polarization p; of 
scattered light for the same conditions as in Fig. 1(a). 


differences as 0—0 in Fig. 1(b) are due to different values of 
ago of particles assumed for the geometric optics case (do 
= 1000 um) and Mie theory (aọ=10 um) calculations. 

The position of the step around 80° [see Fig. 1(b)] depends 
on the particle size and the wavelength. Therefore, white 
light scattered around this angle will be colored. 

The degree of linear polarization p; of scattered light for 
initially unpolarized incident light beam is given in Fig. 2. It 
was calculated using the corresponding equation in Table I. 
The elements ®,, and ®,, [see Eq. (4)] were averaged ac- 
cording to Eq. (12) with the particle size distribution (1). 

We see that p; is rather small both for bubbles and droplets 
in the forward scattering region (below 90°). It increases for 
larger angles and is larger for bubbles everywhere in the 
backward hemisphere except in regions where rainbow 
(around 140°) and glory scattering (around 180°) occurs for 
water droplets. Interestingly, the value of p; for bubbles is 
positive in the backward hemisphere (except at a narrow 
angle range near 180°). It is negative, however, for droplets 
at scattering angles larger than 155°. It becomes positive for 
droplets again around the backward ( 0= 77) scattering direc- 
tions. We see that both bubbles and droplets show the change 
of the sign of p; at the glory scattering region. Generally, the 
angular dependence of p;(0) is much smoother for bubbles 
than for droplets. 

We present the degree of circular polarization p.(0) of 
singly scattered light in Figs. 3(a) and 3(b). It has been cal- 
culated using the corresponding equation in Table I. The el- 
ements ®,, and P44 [see Eqs. (4) and (5)] were averaged 
according to Eq. (12) with the particle size distribution (1). 
The case given in Fig. 3(a) corresponds to incident light that 
is completely right circularly polarized (p.=1). We see that 
both bubbles and droplets do not significantly change the 
degree of polarization for angles smaller than 20°. They both 
convert incident right circularly polarized light (p.=1) into 
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Fig. 3. (a) The same as in Fig. 2, but for the degree of circular polarization 
Pe. (b) The same as in Fig. 2, but for the degree of circular polarization p** . 


left circularly polarized light (p.=—1) at a 180° scattering 
angle. The curve p.(@) is much smoother for bubbles com- 
pared to droplets. It crosses the line p.=0 only once in con- 
trast to droplets, where the degree of circular polarization of 
scattered light p, vanishes at seven scattering angles, there- 
fore, producing partially linearly polarized light (p.=0). 

Bubbles change the sense of rotation of the electric vector 
at angles larger than 105° (p.<0). Droplets are character- 
ized by different signs of p.(@) depending on the specific 
scattering angle interval in the backward hemisphere. 

The degree of circular polarization p*(0) of scattered 
light for incident linearly polarized light (the azimuth angle 
w=45°) is given in Fig. 3(b). It was calculated using the 
corresponding equation in Table I. The elements ®,, and 
@3,4 [see Eqs. (4) and (5)] were averaged as before. 

The value of p* is close to zero for bubbles and it is in the 
range [—0.65,0.4] for water droplets in the backward hemi- 
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sphere. On the other hand, the value p* for droplets is close 
to zero in the forward scattering hemisphere in comparison 
to larger values of p* for bubbles [see Fig. 3(b)]. Larger 
values of p* mean that bubbles can easily convert linearly 
polarized light in the circular polarization mode in the for- 
ward scattering region in contrast to droplets. Droplets give 
large negative values of p* at the rainbow angle. 

The values of the asymmetry parameters g 
=0.25f$p(6)sin20d6@ for both droplets and bubbles are 
close (~0.85). This means that thick layers of both clouds 
and bubbles of the same optical thickness will have similar 
values of the diffuse transmittance and reflection 
coefficients. Such a behavior is in contrast to the differences 
in singly scattered light characteristics shown in Figs. 1-3. 


IV. THE ENTROPY OF SCATTERED LIGHT 


We briefly discuss the difference in entropy of light beams 
scattered by bubbles and droplets. The maximum entropy 
principle is a central organizing principle of statistical phys- 
ics. The entropy of a radiation field is considered in Refs. 7 
and 8. The concept of entropy allows us to introduce a ther- 
modynamic approach to the problems of polarization 
optics.*? In particular, entropy calculations play an important 
role in the analysis of the irreversible evolution of the mul- 
tiply scattered light by a random medium.'° 

The entropy s of the radiation field is defined as!! 


s=—Tr(D InD), (13) 


where Tr means the trace operation and D is the density 
matrix. Equation (13) is a standard definition used in statis- 
tical mechanics. Equation (13) can be reduced to® 


s=—In(a*B*), (14) 


where a=(1+II)/2 and B=(1—II)/2 are eigenvalues® of 
the density matrix. The value of II gives the degree of total 
polarization, defined as VS34+S34+S7/S, .8 Note that Eq. 
(14) also can be written as 


s=—In[0.5(1 +I1)2(1 — I)®]. (15) 


Consider the illumination of a spherical polydispersion 
by a linearly polarized light beam with the azimuth angle 
y= —45°. Then we have (see Table I): 


T= Vp? + p2+p*?. (16) 


IT is exactly equal to unity for monodispersed spheres, which 
can be easily shown using Eqs. (4), (5), and Table I. There- 
fore, the difference, A=II—1, shows the influence of poly- 
dispersity on the value of the degree of total polarization I. 
Note that we have for monodispersed spheres [see Eq. (15)] 
s=0, which means that single scattering by monodispersed 
spheres does not produce entropy of scattering beams. The 
entropy is due to the polydispersity of the medium or mul- 
tiple scattering effects. 

We have for a completely polarized light beam (II=1), 
s=0, independent of the specific type of polarization (linear, 
elliptic or circular). For unpolarized beams (II=0), the en- 
tropy is a maximum and, therefore, s=In 2. 

The function II(@) is shown in Fig. 4. The value of II 
remains larger than 0.9 for bubbles for almost all scattering 
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Fig. 4. The same as in Fig. 2, but for the degree of polarization II. 


angles (except the glory scattering region, where it drops to 
0.4). It is smaller for droplets with three distinctive minima 
at approximately 85°, 118°, and 160°. For the last minimum, 
the value of II drops below 0.1. Therefore, at this angle, 
around 90% of the incident completely polarized light is 
transformed to completely unpolarized light. Note that II 
=1 in the exact forward and backward scattering directions 
because A ,;=A at 0=0 and A;;=—A,, at 0= 7 [see Eqs. 
(4)-(7), (9), and Table I]. 

The angular dependence of the entropy of scattered light is 
presented in Fig. 5, which was derived from data given in 
Fig. 4 using Eq. (15). We note that Fig. 4 corresponds to a 
completely polarized incident light beam (s=0). The en- 
tropy of scattered beam is not equal to zero, however. The 
entropy production As=s,—s,, where sı gives the entropy 
of incident beam and s, is the entropy of scattered beam, 
generally increases with the scattering angle. It drops at the 
rainbow (for droplets) and glory regions due to the high po- 
larizing ability of particles in these regions. The entropy pro- 
duction is equal to zero at exact forward and backward scat- 
tering regions. Spherical particles do not depolarize incident 
completely polarized light in these directions. 


V. CONCLUSION 


In this paper, the author has contrasted the intensity and 
polarization characteristics of scattered light beams for en- 
sembles of droplets in air and air bubbles in water. The scat- 
tering by such systems differs considerably even if the scat- 
ters have exactly the same size. Thus, results obtained for 
droplets cannot be generalized to bubbly media and vice 
versa. The calculations were performed for a single particle 
size distribution f(a). However, our results are valid [at 
least, qualitatively, see Fig. 1(b)] for all nonabsorbing par- 
ticles with a>). They can be used as input for the solution 
of the vector radiative transfer equation for cloudy and bub- 
bly media. 
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Fig. 5. The same as in Fig. 2 but for the scattered light beam entropy s. 


The angular dependence of the entropy of scattered light 
beam for bubbles and water droplets has been studied. Drop- 
lets are far more effective producers of entropy in compari- 
son to bubbles. The use of the notion of entropy and espe- 
cially the principle of its maximum can lead to further 
advances in the field of polarization optics. In particular, it 
may be the goal of optical design to construct a measurement 
system with As=0.° 

The author emphasizes that the Mie theory is not normally 
discussed in the graduate or undergraduate curriculum. How- 
ever, because light scattering by a sphere has such a long list 
of applications,“ physicists should be aware of the subject. In 
particular, light scattering by aerosols and clouds is one of 
the central issues of modern satellite remote sensing of the 
terrestrial troposphere,‘ and is important, for example, for 
global warming issues. The extension of the research given 
here (especially as far as the entropy of the scattered light is 
concerned) would make an excellent student project and help 
students learn about the field of light scattering by small 
macroscopic particles. 
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